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Abstract: 

We review the resonance gas formalism of hadron thermodynamics and recall that an 
exponential increase of the resonance spectrum leads to a limiting temperature of hadronic 
matter. We then show that the number p(n) of ordered partitions of an integer n grows 
exponentially with n and satisfies the integer counterpart of the statistical bootstrap 
equation. Considering the set of all partitions as a Gibbs ensemble provides a partition 
thermodynamics which is also governed by a limiting temperature, determined by the 
combinatorial structure of the problem. Further associating intrinsic quantum numbers 
to integers results in a phase diagram equivalent to that found in QCD for hadronic matter 
as function of temperature and baryochemical potential. 



1. Introduction 

The most important observation in particle physics around 1960 was that the number 
of different species of so-called elementary particles seemed to grow without limit. Hadron- 
hadron collisions produced more and more resonant hadronic states of increasing masses. 
This phenomenon triggered two different theoretical approaches. 

The more conventional idea, classical reductionism as pursued in atomism since an- 
tiquity, proposed that there must be a smaller number of more elementary objects, which 
then bind to form the observed hadrons as composite states. This approach, as we know, 
ultimately led to quantum chromodynamics as the fundamental theory of strong interac- 
tions and thus once more proved most successful. 

A second, truly novel approach asked what such an increase of states would lead 
to in the thermodynamics of strongly interacting matter [Tj. Both the question and the 
answer: the existence of an ultimate temperature of hadronic matter, are the contribution 
of Rolf Hagedorn. We know today that strong interaction thermodynamics leads to critical 
behavior, to a phase transition in which hadronic matter turns into a plasma of deconfined 
quarks and gluons. Statistical QCD has confirmed this and led to a detailed picture of 
this novel phenomenon. Let us see how it comes about and what is its hadronic basis. 

2. Hadron Thermodynamics 

Consider an ideal gas of identical neutral scalar particles of mass m contained in a 
box of volume V, assuming Boltzmann statistics. The grand canonical partition function 
of this system is given by 




(1) 



leading to 
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For temperatures T 3> m , the energy density of the system becomes 
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the particle density 
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and the average energy per particle 



lu ~ 3 T. 



(5) 
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Hence an increase of energy of the system has three consequences: it leads to 

• a higher temperature, 

• more constituents, and 

• more energetic constituents. 

If we now consider an interacting gas of basic hadrons and introduce resonance forma- 
tion as the fundamental feature of their dynamics, we can approximate the interacting 
medium as an ideal gas of the possible resonance species 0. The partition function of 
this resonance gas is 

\nZ(T,V)=J2^p(m l )K 2 (^) (6) 

i 

where the sum begins with the stable ground state m and then includes the possible 
resonances m,,z = 1,2, ... with weights p(rrii) relative to tuq. Clearly the crucial question 
here is how to specify p(rrii), how many states there are of mass ra%. It is only at this 
point that hadron dynamics enters. 

The simplest answer is that p(m) is obtained just from combinatorics, and it immedi- 
ately gives us Hagedorn's statistical bootstrap model. Hagedorn assumed that "fireballs 
consist of fireballs, which consist of fireballs...". In a more modern form, one would as- 
sume that resonance formation and decay follow a self-similar pattern. In any case, the 
defining equation for p(m) is 



^ AT' 

N iv • 



V'o 



p(m,V ) = 5(m—m ) + 

N-l N 

J Y[ \ dm i p( m i) d 3 Pi] 5 4 (S^ -p). (7) 



i=i 



(2tt) 3 

It was solved analytically by W. Nahm giving 

p(m, Vq) = const. m~ 3 exp{m/Tff} 



The density of states thus increases exponentially in m, with a coefficient T H determined 
by 

Y^f(m /T H ) 2 K 2 (m /T H ) = 21n2 - 1, (9) 

in terms of two parameters Vq and m . Hagedorn assumed that the composition volume 
Vq, specifying the intrinsic range of strong interactions, is determined by the inverse pion 
mass as scale, Vq ~ (47r/3)m~ 3 . This leads to a temperature T H ~ 150 MeV. It should be 
emphasized, however, that this is just one possible way to proceed. In the limit m — > 0, 
eq. © gives 

T H = [vr 2 (21n2 - l)] 1 / 3 V ' 1/3 ~ l/r h , (10) 

where Vq = (47r/3)r^ and denotes the range of strong interactions. With ~ 1 fm, 
we thus have Th — 200 MeV. From this it is evident that the Hagedorn temperature 
persists in the chiral limit and is in fact only weakly dependent on m n , provided the 
strong interaction scale V is kept fixed. 



3 



If we now replace the sum in the resonance gas partition function © by an integral 
and insert the exponentially growing mass spectrum (J8J), 



1/7 r m- 
\nZ(T,V) ~ V — dm m 2 p{m t ) K 2 {^) 



we obtain a divergence for all T > T# : in other words, Tjj is the ultimate temperature of 
hadronic matter. In contrast to what we found above, an increase of energy now leads to 

• a fixed temperature limit, T — ► T#, 

• the momenta of the constituents do not continue to increase, and 

• more and more species of ever heavier particles appear. 

We thus obtain a new, non-kinetic way to use energy, increasing the number of species, 
not the momentum per particle. 

Hagedorn originally interpreted Tjj as the highest possible temperature of strongly 
interacting matter. Somewhat later Cabibbo and Parisi [3] pointed out that the resonance 
gas partition function (|11J1 in fact does not diverge at T = T#, while its higher derivatives 
do. Such behavior occurs at phase transition points, and so Tjj is "only" a critical 
temperature of strongly interacting matter. It is clear now that Tjj indeed defines the 
transition from hadronic matter to a quark-gluon plasma. Hadron physics alone can only 
specify its inherent limit; to go beyond this limit, we need QCD. 

The crucial feature leading to the observed limit of hadron physics at T# is the ex- 
ponential increase in the number of hadronic states. In the following section, we want to 
study the origin of such an increase and look somewhat closer at the nature of T#. 



J dm m 3 / 2 exp{— m 



T T, 



H 



3. Partition Thermodynamics 

To arrive at the simplest possible problem leading to an exponentially increasing num- 
ber of 'states', we consider the number p(n) of ordered partitions of an integer n into 
integers. To illustrate: for n = 3, we have the ordered partitions 3, 2+1, 1+2, 1+1+1, so 
that 

p{n = 3) = 4 = 2 n -\ (12) 

Similarly, n = 4 gives the partitions 4, 3+1, 2+2, 2+1+1, 1+3, 1+2+1, 1+1+2, 1+1+1+1, 
leading to 

p(n = 4) = 8 = 2 n -\ (13) 
This solution can be shown to hold in fact for all n (see Appendix), i.e., 

p(n) = 2 11 - 1 = - exp{nln2}. (14) 

We note here for completeness that calculating the number q(n) of unordered partitions 
of an integer n (i.e., not counting permutations) is more difficult and solved only asymp- 
totically 0: 

g(n) = ^eW^}[i+0(i^)] (15) 
4\/3 n 1 n 1 ^ J 
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We will not make use of this result here. 

The problem of ordered partitions can also be solved for restricted or extended cases. 
For example, if we allow only odd integers in the partitions, we get for large n 

p(n)^e nhK (16) 

with k = (1 + V / 5)/2 — 1.62 < 2. On the other hand, if we introduce further degrees of 
freedom, e.g., giving each integer an intrinsic quantum number ('spin') which can take on 
the values ±1, the exact solution is 

pin) = 2n n ~ l = - e {nlnK} , (17) 

now with k = 3 > 2. Hence we have as solution to the generalized ordered partition 
problem the form (see Appendix) 

p(n) oc n n = e {nlnK} (18) 

with k = 2 for the "standard" case and in general oo > k > 1. 

We thus obtain an exponential increase in the number of possible partitions. How is 
this related to the statistical bootstrap condition 0? Consider a density of states defined 
by a bootstrap equation for integers, 

n k 

p{n) = 6{n-l) + E 7J II P( n i) S&m-n). (19) 

k=2 fc - i=l 

Its solution is just the number of partitions of n, 

p{n) = z p(n) (20) 

up to a normalization of order unity (for the standard case k = 2, z ~ 1.25). We can thus 
conclude 

• the statistical bootstrap equation is an integral formulation of a partition problem; 

• the exponential increase in the number of states is of combinatoric origin, and 

• the exponential coefficient k is determined by the combinatoric structure. 

Let us therefore consider k in more detail. In the standard case, see eq. (JH|), the expo- 
nential increase in the number of partitions is determined by In 2; in eq. 0, we had seen 
that in the actual statistical bootstrap, Tjj is also determined by In 2, but in a somewhat 
more complex way due to the presence of momentum degrees of freedom in addition to 
the hadron masses. 

To understand the meaning of k in a thermal context, we construct a statistical me- 
chanics of partitions, considering each partition as a point in the 'phase space' of all 
partitions of fixed n. The set of all partitions thus forms a Gibbs ensemble, and we define 
the partition entropy as 

S(n) = lnp(n) ~ nlnn, (21) 
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the specific partition entropy as 



s(n) 



S(n) 



n 



and the partition temperature 6 as 



1 

6 



dS{n) 
dn 



In k 



In k. 



(22) 



(23) 



The parameter In k determining the increase of the number of partitions of n is thus the 
partition temperature 0, which here for large n coincides with the inverse of the specific 
partition entropy. 

Next we turn to the mentioned generalized form, in which each integer has an 
intrinsic 'spin' variable cr, = ±1, in order to see what the effect of conserved quantum 
numbers is. We define the total spin for a given partition of n as 



a T (n 



E 



(24) 



where the sum runs over all k < n terms of the partitions. We now calculate the number 
p(n, <7t) of partitions of n at fixed or- For oy = n we have p(n, n) = 1, since only 
the partition n = 1 + 1 + ... + 1 of all spins +1 can give a total spin equal to n. For 
a? — we find p(n, 0) ~ 3 n_1 , since total spin zero constitutes the dominant part of all 
configurations. The specific partition entropy s(n, cy) is shown in Fig.^ It was calculated 
numerically for values of n up to 20000 (see Appendix) and it was found to depend only on 
x = <Jt/u. Moreover, our calculations show no appreciable difference between n = 10000 
and n = 20000, so that the latter is in good approximation the asymptotic form of the 
curve. It can be fitted by 

s(x) = [1 -x a } b \n3, (25) 
with a=2.1, b=0.736, and is shown in Fig. Q varying from s(0) = In 3 to s(l) = 0. 



s(n,o T ) 




Figure 1: Specific partition entropy as function of x = ax/n 



The partition temperature now becomes 

1 dS(n,a T ) d[ns(x)] 



9(x) 



dn 



dn 



s(x) 



x- 



ds(x) 
dx 



(26) 
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and is thus no longer equivalent to the inverse specific partition entropy: for x — > 1, both 
s(x) and Q(x) vanish. Through (J26|) we can determine Q(x) from the function s(x) (see 
eq. (125(0 that we have previously calculated. The resulting form is shown in Fig. |21 it 
decreases from = 1/ In 3 at £ = to at £ = 1. 








Figure 2: Partition temperature as function of x = ax/n 

The constituents of our partition thermodynamics are the different integers occurring 
in the partitions; they are the resonances of the resonance gas. It is thus interesting to ask 
for the relative abundance R n>aT (a, b) of two integers a, b < n in the set of all partitions 
of n. In the appendix, we calculate the total occurrence h(n, a, ax) of an integer a in the 
partitions of n at fixed <jt- Given this, the desired relative abundance is defined as 

h(n, a, cry) 



R 



n.aT 



[a,b) 



h(n, b, ax) 



(27) 



For a, b <C n it turns out that for large n, i? n o - T (a, b) depends only on the ratio ax/n and 
has the form 



Rn,<T T {a,b) ~ e 



-(a-b)/Q{a T /n) 



(2? 



where G is precisely the partition temperature shown in Fig. |21 We thus see that the 
resonance gas feature of species abundances specified by one universal temperature in fact 
has its origin in the underlying combinatoric structure. 

So far, we have considered partition thermodynamics in what might be called a grand 
microcanonical formulation: the number of terms in the different partitions varies, but 
their sum n and the overall spin ot are specified exactly. To go to a grand canonical 
formulation, we fix n and ot only on the average. The starting point is the grand canonical 
partition function 



/*oo rn 

Z(T 1 fi)= dnexp{—n/T} <io"rexp{— /xctt/T} p(n, ut) 

Jl J-n 



(29) 



with the temperature T and the chemical potential \x the Lagrangian multipliers for n 
and ot- Since p(n, <tt) = exp{ns(crr/ra)}, where s^ar/n) is the specific entropy calculated 
above, we get 



/oo rn 
dn J da T exp{— n 



1 fia T ( i \ 

7^ — 7^ slaT/n) 

T T n y ' ' 



}■ 



(30) 
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The integral (j3*Uj) is not defined for all pairs of values (T, /x): for fixed /x, it exists only for 
those values of the temperature for which 

- - s{a T n) > => T < \/a T , n. 31 

i in 

We note that the function on the right of the last inequality depends only on the variable 
x = ax/n, and not on a? and n separately. Since the condition (|31|) must be valid for any 
possible value of the ratio (7t/ti, with \<jj</ri\ < 1, T has to be smaller than the absolute 
minimum of the function 

f(x^) = ^f, (32) 
s{x) 

in the interval —1 < x < 1. Hence for any /x there exists a limiting value T c (/x) of 
the temperature beyond which the grandcanonical partition function is not defined. In 
particular, for /x = 0, we get T c = 1/ln 3, which is the temperature determined in the 
previous section for the partition problem with total spin ot = 0. We further note that 
|/x| < 1, because otherwise the function f(x,fi) would become negative for some x and 
the inequality (fBTj) could never be satisfied by positive values of the temperature T. By 
varying /x we obtain a curve in the T — /x plane, which defines the existence region for our 
system. Determining numerically the minimum of the f(x, /i) as function of /x, we obtain 
the form shown in Fig. El 

The resulting 'phase diagram' of partition thermodynamics, i.e., the boundary for the 
existence of such thermodynamic systems, is evidently very similar to that expected in 
statistical QCD for the phase boundary of hadronic matter. It thus seems that quantum 
chromodynamics leads to a resonance pattern which is basically of combinatoric nature. 
Moreover, the conservation of baryon number also appears to follow a simple combinatoric 
form resulting from the addition of the baryonic degrees in the resonance composition. 
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Figure 3: Phase boundary for partition thermodynamics 
4. Conclusions 

We have shown that the limit of hadron physics obtained in resonance gas thermo- 
dynamics with an exponentially increasing mass spectrum is due to the combinatorics of 
resonance formation. In particular, the simplest partition systems provide 
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• a limiting temperature determined solely by the combinatoric structure of the sys- 
tem; 

• a reduction of this limiting temperature for conserved intrinsic quantum numbers 
and hence a phase boundary for partition thermodynamics; 

• a relative abundance of constituents governed by the limiting temperature. 

Hence, if all combinations of n hadrons are resonances and the degeneracy of each reso- 
nance is determined by the number of composition patterns, an energy input eventually 
goes into species formation, not into faster constituents. Comparing such a resonance gas 
to a conventional pion gas, we thus find for the dependence on energy density e 

pion gas resonance gas 

~ e 3/4 n res ~ e 

LU n ~ e 1//4 uj Tes ~ const. 

Here n denotes the density of constituents, ui the average energy of a constituent. 

We close with a rather speculative question, concerning the relation of our considera- 
tions to nuclear collisions. In the experimental study of such interactions, it is found that 
the relative abundances of hadron species are essentially determined by two parameters, 
the temperature T of a resonance gas and the baryochemical potential /i specifying its 
overall baryon density jHj. For high energy collisions at /i = 0, T ~ 175 MeV; with in- 
creasing /i, T(n) decreases in a form quite similar to that shown in Fig. El What can we 
learn from this? There seem to be two possible scenarios: 

• Each collision produces a thermal system and thus corresponds to an ensemble of 
many partitions: nuclear collisions produce matter. 

• Each collision is one partition, and only the ensemble over many collisions forms a 
thermal ensemble: nuclear collisions simulate matter. 

Which of these two is realized in nature appears to be one of the most crucial questions 
for the experimental high energy heavy ion program. 
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Appendix 

A. Standard Case 

To determine the number of ordered partitions of an integer n into integers ki < n, we 
start from the number of partitions containing m terms. This is equivalent to the number 
of divisions of a segment of length n units into m parts. Each partition of this type is 
constructed by choosing m — 1 of the first n — 1 integers as end points for the m segments 
dividing [0, n]. This can be done in 

p(n, m) = I n 1 ) (33) 
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possible ways. By summing over all possible values of m, we obtain 

n n ( n - \\ 

p(n)=EM = E _i = 2 ^ 

m=l m=l \ m V 

Next, we determine the abundance of a given integer a in the set of all partitions. Suppose 
we have a partition of m terms with at least one occurrence of the number a < n. By 
changing the position of a within the partition, we get m possible partitions containing the 
same number of terms, with order of the m — 1 other terms unchanged. Hence the number 
of occurrences of a in partitions with m summands is just m times the number of partitions 
of n— a into m— 1 terms. Among the partitions of n—a into m— 1 terms there will in general 
be other a's, and in this case each partition gives a contribution larger than one to the 
occurrence of a. However, this larger contribution is exactly compensated by a reduction 
of the number of corresponding partitions, since the different a are indistinguishable and 
can be freely permuted without generating new partitions. To get the total occurrence 
h(n, a) of a in all partitions of n we sum over all possible values of m and obtain 

( n -a-l\ 

{U,a)= £ m { m-2 J' (35) 

Defining I = m — 2, we get 



n-a-l /„ . lN n-a-l /„ q — \\ n ~<^ 1 ( n - CL 



(z+ 2) r ; >2E : +EM 7 • (36) 



and from this 



h{n, a) = {n-a-l) T^ 2 + 2""" = — - T' a (37) 



for the abundance of a in the partitions of n. 

B. Intrinsic Degrees of Freedom 

We now assume that each integer carries an intrinsic degree of freedom, an integer 
'spin' a, giving it 2a degenerate states; in the example above, we had a = 1 and hence 
the two states ±1. It is obvious that such a spin leads to an increase in the number of 
partitions; each partition is now labeled by n and by the total spin a? = J2k a k of its 
summands. If we do not impose any constraint on the total spin, it is trivial to calculate 
the total number of partitions. There are now 2a copies of each number, so that in a 
partition with m terms, each summand contributes an additional factor 2a. Eq. (}3*3*|) thus 
becomes 

p a (n,m) = (2a) m ( n ~]) (38) 



,m 



so that 



Pa{n)= ^j9 CT (n,m)= £ (2a)™ ( _ 1 ) = (2a) (2a + (39) 
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The problem is more interesting if ax has a well defined value, which is equivalent to 
requiring the conservation of a spin, isospin or baryon number for a system of particles. In 
this case, it is not possible to compute the total number of partitions and the abundances 
in closed form. However, at least for the special case a = 1, we can derive formulas 
that allow a numerical computation of the quantities of interest. Since the problem is 
completely symmetric under a change of sign of the spin, we consider positive values of 
the total spin. 

We start by asking how many partitions of m terms there are with a total spin ot- 
The summands of the partitions carry spin +1 or —1. If there are k numbers with spin 
+1 there must then be m — k numbers with spin —1. In order to have a total spin o>r, 
there is only one possibility: (m + a?)/ 2 terms must carry spin +1 and (m — oy)/2 terms 
spin —1. Hence the required number of partitions at fixed m is 

p(n,m,a T )= (. ™ Wo ] T ~\ ) (40) 



(m + a T )/2 
and the total number becomes 

P[n '° T) = „£/ ( WT) = L T ((m Zr)ft) (m- I} (41) 

Note that the sum over m is constrained by the parity of the total spin: if ot is even 
(odd), m must be even (odd). Hence the sums run over those values of m between ot 
and n which have the same parity as <jt- 

To obtain the abundances, we only have to insert in Eq. ()35j) the factor due to the 
spin degeneracy of the partitions, given in Eq. (jUJ), leading to 

h(n,a,a T )= £ ml, , _ Wo _ „ . (42) 



(m + a T )/2 



Eqs. (j41j) and (|4*2^) can easily be evaluated using Mathematica. 
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